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Abstract.  We  prove  the  existence  of  a  weak  solution  to  a  two-dimensional  resonant  3x3  system  of 
conservation  laws  with  BV  initial  data.  Due  to  possible  resonance  (coinciding  eigenvalues),  spatial  BV 
estimates  are  in  general  not  available.  Instead,  we  use  an  entropy  dissipation  bound  combined  with 
the  time  translation  invariance  property  of  the  system  to  prove  existence  based  on  a  two-dimensional 
compensated  compactness  argument  adapted  from  [36].  Existence  is  proved  under  the  assumption  that 
the  flux  functions  in  the  two  directions  are  linearly  independent. 


1.  Introduction 


This  paper  studies  certain  two-dimensional  resonant  3x3  systems  of  conservation  laws  of  the  form 

kt  =  0,  h  =  0, 

Ut  fi^kj  u)x  11.)^  —  0, 

which  are  augmented  with  L°°  n  BV  initial  data 


(1.2)  k\t^Q  =  k{x,y),  l\t^Q  =  l{x,y),  u\t=o  =  uo{x ,  y) . 

The  goal  is  to  prove  that  there  exists  a  weak  solution  to  (1.1)-(1.2). 

In  recent  years  the  one-dimensional  version  of  the  above  system, 

kt  =  0, 

ut  +  f{k,u)x  =  G, 


has  received  a  considerable  amount  of  attention.  This  system  may  be  viewed  as  an  alternative  way  of 
writing  a  scalar  conservation  law  with  a  discontinuous  flux,  namely 


(1.4) 


ut  +  f{k{x),u)x  =  0. 


Equations  like  (1.4)  occur  in  a  variety  of  applications,  including  flow  in  porous  media,  sedimentation 
processes,  traffic  flow,  radar  shape-from-shading  problems,  blood  flow,  and  gas  flow  in  a  variable  duct. 

If  k{x)  is  a  smooth  function,  Kruzkov’s  theory  [21]  tells  us  that  there  exists  a  unique  entropy  solution 
to  the  initial  value  problem  for  (1.4),  for  general  flux  functions  /.  The  scalar  Kruzkov  theory  does  not 
apply  when  k{x)  is  discontinuous.  Instead  it  proves  useful  to  rewrite  (1.4)  as  a  2  x  2  system  of  equations 
(1.3),  which  makes  it  possible  to  apply  ideas  from  the  theory  of  systems  of  conservation  laws. 

As  a  starting  point,  it  is  necessary  to  introduce  conditions  on  the  flux  /(fc,  u)  that  guarantee  that 
solutions  stay  uniformly  bounded.  For  example,  one  can  require  /(fc,  a)  =  f{k,  6)  =  0  for  all  k,  which  in 
fact  implies  that  the  interval  [a,  6]  C  K  becomes  an  invariant  region.  The  system  (1.4)  has  two  eigenvalues, 
namely  Ai  =  0  and  A2  =  fu{k,u).  Consequently,  if  fu{k,u)  vanishes  for  some  value  of  {k,u),  then  (1.4) 
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is  nonstrictly  hyperbolic  and  experiences  so-called  nonlinear  resonant  behavior,  which  implies  that  wave 
interactions  are  more  complicated  than  in  strictly  hyperbolic  systems.  As  a  matter  of  fact,  one  cannot 
expect  to  bound  the  total  variation  of  the  conserved  quantities  directly,  but  only  when  measured  under 
a  certain  singular  mapping.  A  singular  mapping  that  is  relevant  for  (1.3)  is 


'i>{k,u)  =  l  \fu{k,0\dt 


If  {u^}pyQ  is  a  sequence  of  ” reasonable”  approximate  solutions  of  (1.3),  then  one  proves  that  the  total 
variation  of  the  transformed  quantity  :=  'I'(fc,  m^)  is  bounded  independently  of  p.  Helly’s  theorem 
then  gives  convergence  (along  a  subsequence)  of  as  p  i  0.  Since  the  continuous  mapping  u  ’i’{k,  u) 
is  one-to-one,  also  converges. 

A  singular  mapping  was  used  first  by  Temple  [39]  to  establish  convergence  of  the  Glimm  scheme  (and 
thereby  the  existence  of  a  weak  solution)  for  a  2  x  2  resonant  system  of  conservation  laws  modeling 
the  displacement  of  oil  in  a  reservoir  by  water  and  polymer,  which  is  now  known  to  be  equivalent  to  a 
conservation  law  with  a  discontinuous  coefficient  (see,  e.g.,  [20]).  Since  then  the  singular  mapping  ap¬ 
proach  has  been  used  and  adapted  by  great  many  authors  to  prove  existence  of  weak  solutions  to  resonant 
systems  of  conservation  laws/scalar  conservation  laws  with  discontinuous  flux  functions,  by  establishing 
convergence  of  various  approximations  schemes  (Glimm  and  Godunov  schemes,  front  tracking,  upwind 
and  central  type  schemes,  vanishing  viscosity/smoothing  method,  ...),  see  (the  list  is  far  from  being 
complete)  [1,  2,  3,  4,  5,  6,  8,  9,  10,  11,  12,  18,  19,  20,  22,  23,  29,  35,  40,  41].  Similar  ideas  have  been  used 
also  in  the  context  of  degenerate  parabolic  equations  [15].  Regarding  uniqueness  and  entropy  conditions 
for  scalar  conservation  laws  with  discontinuous  coefficients,  see  [16,  17]  and  the  references  therein. 

As  an  alternative  to  the  singular  mapping  approach,  the  papers  [13,  14,  17]  has  suggested  to  use  the 
compensated  compactness  method  and  ’’scalar  entropies”  for  the  convergence  analysis  of  approximate 
solutions.  The  results  obtained  with  this  approach  are  more  general  (and  to  some  extent  the  proofs  are 
easier)  than  those  obtained  with  the  singular  mapping  approach. 

All  the  papers  up  to  now  have  addressed  the  one-dimensional  case.  The  aim  of  the  present  paper  is 
to  take  a  first  look  at  the  multi-dimensional  case,  which  is  completely  unexplored.  More  precisely,  we 
will  prove  the  existence  of  at  least  one  weak  solution  to  the  initial  value  problem  for  the  two-dimensional 
system  (1.1). 

Our  existence  proof  is  based  on  studying  the  ”(£,<5)  i  (O,!*)  limit”  of  classical  solutions  of  the 
uniformly  parabolic  equation 


e,(5 


+  f{k  ,  )x  +  9{l  7  )y  —  ^  {'^xx  +  ) 


£  >  0,  5  >  0, 


where  k^ ,  converge  to  k,  I  in  Ljqj,(R^),  respectively,  as  5  J,  0. 

Observe  that  we  are  essentially  considering  a  scalar  approximation  scheme  for  (1.1),  see  [5,  6,  15,  13, 
14,  17,  40,  41]  for  other  scalar  approximation  schemes  for  one-dimensional  discontinuous  flux  problems. 

Although  spatial  BV  bounds  are  out  of  reach,  we  still  have  a  time  translation  invariance  property 
at  our  disposal,  which,  together  with  the  assumption  of  BV  initial  data,  implies  that  is  uniformly 
bounded  in  L^.  Gonsider  three  functions  F{k,u),  G{l,u),  H{k,l,u)  defined  by 


Fu  =  Gu  =  {Puf,  Hu  =  fu9u 


We  prove,  at  least  under  the  assumption  that  e  and  5  are  of  comparable  size,  that  the  two  sequences 

F{k{x,  y),  +  H{k{x,  y),  l{x,  y),u^’^)y 


and 

H{k{x,  y),  l{x,  y),u^'^)x  +  G{k{x,  y),u^'^)y 
are  compact  in  ’^(R^),  for  each  fixed  t  >  0. 
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The  crux  of  the  convergence  analysis  is  then  to  prove  that  the  above  compactness  is  suffi¬ 

cient  to  establish  a  ’’two-dimensional”  compensated  compactness  argument  in  the  spirit  of  the  classical 
Tartar-Murat  results  for  one-dimensional  conservation  laws,  [26,  27,  28,  37,  38].  Here  we  follow  the 
recent  two-dimensional  compensated  compactness  framework  developed  in  Tadmor  et.  al.  [36]  for  non¬ 
linear  conservation  laws.  We  extend  their  results  to  the  case  involving  additional  discontinuous  ’’variable 
coefficients”.  Accordingly,  we  make  the  nonlinearity  assumption  that  for  each  fixed  k,l  the  functions 
u  1-^-  fu{k,u)  and  u  i-^-  are  almost  everywhere  linearly  independent  (see  (2.4)  in  the  next  section 

for  a  precise  statement) .  Our  main  existence  result  is  based  on  an  application  of  the  two-dimensional  com¬ 
pensated  compactness  lemma  with  ’’variable  coefficients”  —  lemma  3.2  stated  in  Section  3  below.  Granted 
the  nonlinearity  assumption,  it  then  yields  that  (a  subsequence  of)  u'^’^{-,-,t)  converges  in  Ljqj,(R^)  to 
a  bounded  function  u{-,-,t),  for  a.e.  t  >  0.  Since  is  uniformly  Lipschitz  continuous  in  time  we 
obtain,  in  Section  4  below,  our  main  Theorem  2.1,  stating  that  ^  u  in  x  K_|_)  and  that  the 

limit  function  m  is  a  weak  solution  of  (1.1)-(1.2). 

Although  we  have  chosen  to  analyze  the  vanishing  viscosity/smoothing  method,  the  techniques  used 
here  for  that  purpose  can  also  be  applied  to  various  numerical  schemes,  including  appropriate  two- 
dimensional  versions  of  the  scalar  finite  difference  schemes  studied  in  [13,  17,  40,  41]. 


2.  Assumptions  and  statement  of  main  results 

We  start  by  listing  the  assumptions  on  the  initial  conditions  Ug  and  the  fluxes  k,  I,  /,  g  that  are  needed 
for  the  existence  result. 

Regarding  the  initial  function  we  assume 

(2.1)  Ug  G  L°°(R^)  n  BV{R'^),  a  <  Ug  <  b  for  a.e.  in 


(2.2) 


For  the  discontinuous  coefficients  k,l  :R.‘^  we  assume 

( k  I  G  ^ 


)n  w(k2), 

a  <  k,l  <  P  a.e.  in 

For  the  flux  functions  /,  g  :  [a,  /3]  x  [a,  b]  - 


we  assume 


(2.3) 


u  1-^  f{k,  u),  u  1-^  g{l,  u)  G  C^[a,  b]  for  all  fc,  I  G  [a,  /3]; 
k  1-^  f{k,  u),l  1-^  g{l,  u)  G  [a,  /3]  for  all  u  G  [a,  6]. 


Moreover,  we  make  the  nonlinearity  assumption  which  excludes  the  possibility  of  5i/(fc,  u)  ^2g{l,  u) 
being  an  affine  function  (in  u)  on  any  nontrivial  interval  for  all  k,  I  G  [a,  /3], 

(2.4a)  V  1^1  =  1  and  k,  I  G  [a,  /3]  :  ^if{k,  •)  -I-  ^2g{l,  •)  ^  affine  function  on  any  nontrivial  interval. 


In  its  slightly  stronger  version,  this  assumption  requires  that  fu{k,-)  and  gu{l,-)  are  a.e.  linearly  inde¬ 
pendent  so  that  the  symbol  s{^,k,l,u)  ■=  ^ifu{k,u)  ^2gu{l,u)  satisfies 

(2.4b)  V|^|  =  l  and  k,lG[a,P]:  meas{u  \  s{^,  k,l,u)  =  0}  =  0. 

This  is  a  straightforward  generalization  of  the  notion  of  nonlinearity  found  in  [24] ,  in  their  study  of  kinetic 
formulations  for  nonlinear  conservation  laws. 

Finally,  we  need  to  know  that  our  approximate  solutions  stay  uniformly  bounded.  For  example,  this 
is  ensured  by  the  assumption 

(2.5)  f{k,  a),  f{k,  b),  g{l,  a),  g{l,  6)  =  0  for  all  k,  I  G  [a,  /3], 

which  implies  that  the  interval  [a,  6]  becomes  an  invariant  region. 

We  are  now  ready  to  state  our  main  result. 
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Theorem  2.1.  Suppose  (2.1),  (2.2),  (2.3),  (2.4),  and  (2.5)  hold.  Then  there  exists  a  weak  solution 

u  G  L“(R2  X  R+)  n  Up{R+;L\m.‘^)) 


of  the  initial  value  problem  (1.1)-(1.2),  that  is,  u  satisfies 


(2.6) 


/  {u4>t  +  f{k{x,y),u)4>,,;  +  g{l{x,y),u)4>y)  dxdt 
+  /  uo{x)4>{x,Q)  dx  =  Q,  V(/)  G  I?(K^  X  [0,  oo)). 

Jr 


The  weak  solution  u  can  he  constructed  as  a  strong  Lj(,j.(R^  x  ]R_|_)  limit  of  classical  solutions  u'^’^  of 
uniformly  parabolic  problems  with  smooth  coefficients.  Moreover,  if  v  is  a  another  weak  solution,  with 
initial  data  Vq,  constructed  in  this  way,  then  for  any  t  >  0 


(2-7)  ||w(-,t)  '^’(•5  i)||^i(-R2)  <  ||mo  'i'o|lii(R2) , 

and  thus  the  constructed  limit  u  is  unique. 


The  proof  of  this  theorem  is  given  in  the  following  two  sections.  We  close  this  section  with  the  following 
summary. 

Remark  2.2.  Stated  differently,  Theorem  2.1  shows  that  there  exists  a  weak  solution  to  the  following 
two-dimensional  scalar  conservation  law  with  discontinuous  coefficients  k,  I  G  L°°(R^)  n  i3y(R^): 

ut  +  f{k{x,  y),u)a;  +  g{l{x,  y),u)y  =  0, 

u\t^Q  =  Wo  G  T“(R2)  n  W(R2). 

We  note  in  passing  that  the  solution  operator  in  this  case  of  discontinuous  ’’variable  coefficients”  is  not 
translation  invariant  in  space  and  hence,  the  L^-contraction  statement  in  (2.7)  does  not  imply  spatial 
BV  compactness. 

Moreover,  if  we  let  St  :  L°°(R)  ^  L°°(R)  denote  the  solution  operator  of  (1.1)-(1.2),  so  that  uft,-)  = 
SfUof)  is  the  weak  solution  constructed  in  Lemma  2.1,  then  by  adapting  standard  arguments,  we  can 
prove  that  the  solution  operator  St  is  compact  with  respect  to  the  norm. 


3.  A  COMPENSATED  COMPACTNESS  LEMMA 

In  this  section  we  prove  a  ’’two-dimensional”  compensated  compactness  lemma.  We  refer  [7,  25,  26, 
27,  28,  37,  38]  for  background  information  on  the  compensated  compactness  theory.  We  start  by  recalling 
the  celebrated  div-curl  lemma. 

Lemma  3.1  (div-curl).  Let  C  R^  6e  an  open  domain  and  let  p  >  0  denote  a  parameter  taking  its 
values  in  a  sequence  which  tends  to  zero.  Suppose  D,  ^  E  in  (L^(n))  and  {divU^j^^g, 

{curlA^I^^P  belong  to  a  compact  subset  ofW^^'^ifil).  Then,  after  extracting  a  subsequence  if  necessary, 
we  have  Dp  ■  EP  ^  D  ■  E  in  V  {LI)  os  p  J,  0. 

The  compensated  compactness  lemma  below  is  tailored  for  two-dimensional  equations,  whose  spatial 
part  involve  discontinuous  coefficients: 

f{k{x,  y),v{x,  y)),,  +  g{l{x,  y),v{x,  y))y. 

If  g{l,  u)  =  g{u)  and  f{k,  u)  =  f{u),  then  the  lemma  below  coincides  with  the  two-dimensional  result 
of  [36,  Theorem  3.1].  If  we  set  p  =  0  then  the  result  coincides  with  Tartar’s  compensated  compactness 
lemma  for  the  one-dimensional  scalar  conservation  with  genuinely  nonlinear  flux  /. 
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Lemma  3.2  (Compensated  compactness).  Let  C  &e  an  open  domain.  Let  k,l,f,g  be  functions 
satisfying  (2.2),  (2.3),  (2.4),  and  (2.5).  Suppose  {v'^ {x ,  y)} is  a  sequence  of  measurable  functions  on 
LI  that  satisfies  the  following  two  conditions: 

1.  There  exist  two  finite  constants  a  <  b  independent  of  p  such  that 

a  <  v^{x,  y)  <b  for  a.e.  {x,  y)  G  LI. 

2.  Let  the  functions  F,  G,  H  be  defined  by 

Fu{k,u)  =  {fu{k,u)f  ,  Gu{l,u)  =  {gu{l,u))'^  ,  Hu{k,l,u)  =  fu{k,u)gu{l,u). 


We  assume  that  the  two  sequences 


(3.1) 


^F  {k{x,  y),  yP)^  +  H  {k{x,  y),  l{x, 

{k{x,y),l{x,y),vP)^  +  G  {l{x,y),vP)y'^ 


belong  to  a  compact  subset  ofW^^^''^(R  x  M+). 

Then,  there  exists  a  subsequence  of  {yP {x ,  y)} that  converges  a.e.  to  a  function  v  G  L°°(R^),  and 
a  <  v{x,  y)  <  b  for  a.e.  {x,  y)  G  LI. 


Proof.  To  simplify  the  notation  let 

FP:=F{k{x,y),vP),  GP  :=  G  {l{x ,  y) ,  yP) ,  HP  :=  H  {k{x,y),l{x,y),vP) , 
and  denote  their  L°°{R  x  ]R_|_)  weak-x  limits  by  F,  G,  H,  respectively.  Introduce  the  vector  fields 

DP  =  (^FP,  HP^  EP  =  (^-GP,  HP^ , 

and  denote  their  respective  L°°{R  x  ]R_|.)  weak-x  limits  by  D,E. 

Thanks  to  (3.1),  we  can  apply  the  div-curl  lemma  to  the  sequences  produce 

DP  ■  EP  =  D  ■  E  a.e.  in  Lt, 


that  is, 

{Hpf  -  FPGP  =  (uf  -FG, 

which  implies 

(3.2)  (^HP -TT)^  -  {FP -F)  {GP -G)  =  Q. 

Fix  c  =  c{x,  y)  G  L°°{Lt).  Following  [36],  we  now  consider  the  function  /  :  [a,  6]  ^  K  defined  by 

I{v)  =  (h  {k{x,  y),  l{x,  y),v)  -  H  {k{x,  y),  l{x,  y),  c)  ) 

-  (^F  {k{x,  y),v)-  F  {k{x,  y),c)^  ■  (g  {l{x,  y),v)  -  G  {l{x,  y),  c)  ^ 

Note  that 

I{vP)  =  (^[HP -H]  +  [H-H{k{x,y),l{x,y),c)]y 

-  ([F^-F]  +  [F-F{k{x,y),c)])  ■  (^[GP-G]  +  [G  -  G{l{x,y),c)]). 
Using  this  and  (3.2),  we  compute 

I{vP)  =  [h  -  H  {k{x,y),l{x,y),c)^  -  {f  -  F  {k{x,y),c)y{G  -  G  (l(a;,  y),  c)^ 


(3.3) 
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By  the  Cauchy-Schwarz  inequality  we  have  for  any  u  S  [a,  6] 

{h  {k{x,  y),  l{x,  y),u)-  H  {k{x,  y),  l{x,  y),  c)^ 

<  (^F  {k{x,  y),u)-  F  {k{x,  y),c)^  (g  {l{x,  y),u)  -  G  {l{x,  y),  c)^ , 

and  hence  /(•)  <  0  with  1(c)  =  0.  Thanks  to  (2.4),  the  Cauchy-Schwarz  inequality  in  (3.4)  is  in  fact  a 
strict  inequality.  This  shows  that  the  function  I(v)  has  a  strict  global  maximum  at  v  =  c. 

Since  F  (k(x,  y),  u)  is  strictly  increasing,  we  can  choose  c  as 

c(x,  y)  :=  {k(x,  y),  F(x,  y))  , 

so  that  (3.3)  becomes 

(3.5)  I(vP)  =  (^H  -  H  (k(x,y),l(x,y),c)^  . 

Since  /(•)  <  0,  we  conclude  that  FI  =  F[  (k(x,y),l(x,y),c),  and  thus  I(vP)  =  0.  In  fact,  we  have 
I{vP)  0  a.e.  in  fl. 

Using  the  fact  that  c  is  a  strict  maximizer  of  I(v),  we  have 

I(x)  <  —Ga,  whenever  |'u  —  c|  >  a, 
for  some  constant  Ga  >  0  that  depends  on  a.  Consequently, 

(3.6)  meas  —  m|  >  a}  <  ^  ff  I{v^{x,y))  dxdy  ^  0  as  p  J,  0. 

J  J  Q,r\\vP—u\>OL 

Since  a  >  0  was  arbitrary,  yP  ^  c  in  measure,  which  in  turn  implies  that  a  subsequence  of  {r'^},5>g 
converges  to  c  a.e.  in  U.  □ 

We  remark  that  the  idea  of  using  the  Cauchy-Schwarz  inequality  along  the  lines  of  (3.4)  for  proving 
strong  compactness  can  be  traced  back  to  [33,  34] 


4.  Proof  of  Theorem  2.1 


/^OO 

^0 


Let  k^ ,  ,  Uq  be  smooth  functions  converging  strongly  to  k,  I,  uq  respectively.  More  precisely,  let  cos  € 

be  a  nonnegative  function  satisfying 


uj(x)  =  0  for  jccj  >  1, 


uj(x)  dx  =  1. 


For  (5  >  0,  let  ujs{x)  =  (|)  and  introduce  the  mollified  functions 

k^  =  U!s*k,  f=ujs-kl. 

We  approximate  the  initial  data  Ug  by  cut-off  and  mollification  as  follows: 

Mg  =  UJS  *(uox^i), 

where  xs(x)  =  1  for  jccj  <  1/6  and  0  otherwise.  In  particular,  we  have  the  estimate 

1 


(4.1) 


K),,  +  K), 


< 

Li(R2)  6 


{ut)y)  dxdy<^\uo\ 


By(R2)  ■ 


Observe  that  for  =  k^,  Ug  and  h  =  k,l,  ug,  we  have  €  G°°(R^)  and 

^  h  a.e.  in  and  in  L^(R^)  for  any  p  G  [1,  oo)  as  <5  J,  0. 
Additionally,  Ug  is  compactly  supported. 
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Now  let  be  the  unique  classical  solution  of  the  uniformly  parabolic  equation 
(4.2)  +  f{k\  +  g{f,  =  £  (<’i  +  ul’^) , 

with  initial  data  =  Uq.  The  proof  proceeds  through  a  series  of  lemmas,  which  in  the  end  show 

that  for  each  t  G  [0,  T]  a  subsequence  of  ■,  t)  converges  a.e.  as  e,  <5  J,  0. 

Our  first  lemma  confirms  the  uniform  bound. 

Lemma  4.1  {L°°  bound).  There  exists  a  constant  C  >  0,  independent  of  £,6,  such  that 


L“(R2) 


<  C,  for  all  t  G  (0,  T). 


Proof.  The  proof  is  standard  and  exploits  assumption  (2.5)  to  conclude  that  a  <  u'^’^{x,y,t)  <  b  for 
a.e.  {x,  y,  t)  G  x  K_|_ .  □ 

Using  that  (4.2)  is  translation  invariant  in  time,  we  can  prove  that  ul’^  is  uniformly  bounded  in 
L“(R+;L1(K2)). 

Lemma  4.2  (Lipschitz  regularity  in  time).  Suppose  the  two  smoothing  parameters  e  and  5  are  kept  in 
balance  in  the  sense  that 

(4.3)  6  =  Ce,  for  some  constant  C  >  0. 

There  is  a  constant  Cq  (which  is  possibly  dependent  on  uq  but  otherwise  is  independent  of  £,5),  such  that 
for  any  t  >  0 


[[  \dtu‘^^^{-,-,t)\  dxdy  <  Cq. 

J  J 


Proof.  To  prove  this,  set  =  ul’^ .  Then  w  satisfies 

+  {fu{k\  +  {gu{l\  =  e  [wli  +  wl'y) . 

Multiplying  by  sign(r(;®’'^)  gives,  in  the  sense  of  distributions, 

+  ifuik^u^'^)  |u;"’^|),  +  {gu{l\u^'^) 


-  csign'(u;^'^)  +  (u;f  )^)  , 


since  fu{k^ ,  u'^’^)w'^’^sign  =  gu{k^ ,u'^’^)w'^’^sign{w'^’^)wy’^  =  0.  Hence 

^ff  \w^’^\{x,y,-)dxdy<0, 
at  J  J  R2 

which,  due  to  (4.1)  and  (4.3),  concludes  the  proof.  □ 

Thanks  to  the  previous  lemma,  we  also  have  uniform  L°°(R+;  L^(R^))  control  over  \/£u%’^  and  ^/eUy’^. 
Lemma  4.3  (Entropy  dissipation  bound).  There  is  a  constant  C,  independent  of  £,6,  such  that 

£  [[  +  {ul'\-,-,t))^^  dxdy  <  C,  for  any  t  >  0. 
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Proof.  Multiplying  (4.2)  by  and  then  integrating  yield 
J j  ^  •,  t))^  +  dx  dy 

=  —  [ [  dx  dy 


fuik^Odfj  /fe(fc^C)d?jfcf(x,2/) 

J  9u{l\0df\  -  (j  9i{l^,0df\ly  dxdy. 


In  view  of  Lemmas  4.1  and  4.2  and  the  BV  regularity  of  the  coefficients,  we  derive  easily  the  uniform 
bound 

-  ^  (lh‘’1Loo(R^.ii(R2))  +  l*lw(R=)  +  l^lw(R=))  >  ^ 


lL°°(R+;Li(*^)) 


for  some  constant  C  that  is  dependent  on  ||u®’'^||  j^oo(][j2)  but  otherwise  is  independent  of  £,5. 


Lemma  4.4  (Pre-compactness  at  each  time  instant).  Suppose  the  two  parameters  e  and  S  are  kept  in 
balance  in  the  sense  that  (4.3)  holds.  With  F ,  G,  and  H  defined  in  Lemma  3.2,  the  two  sequences 

|f  {k{x,  2/),  -f  H  {k{x,  y),  l{x,  y), 

|iL  {k{x,  y),  l{x,  y),  -h  G  {l{x,  y), 

then  belong  to  a  compact  subset  o/ ’^(R^),  for  each  fixed  t  >  0. 

Proof.  Let  4>  =  4>{x,  y)  G  I?(R^),  and,  for  each  fixed  t  >  0,  introduce  the  distribution 

=  JJ  {k{x,y),u^’^)  (j)a;  +  H  {k{x,y),l{x,y)u^’^)  (fy'j  dxdy. 

Let  us  first  write  where 

=  JJ  Jf  {k{x,y),u‘^'^)  -  F  {k^{x,y),u‘'^^)'j(j)a;dxdy 

+  JJ  Jh  {k{x,y),l{x,y),u'^’^)  -  H  {k^{x,y),l\x,y),u^’^)^4>ydxdy, 

4’)  =  JJ  Jf  {k\x,  y),u^'^)  4x  +  H  {k\x,  y),  l\x,  y),u^'^)  (fy)  dx  dy. 

In  what  follows,  we  let  fl  denote  an  arbitrary  but  fixed  bounded  open  subset  of  M^.  With  (f  G  WQ^’^(n), 
we  have  by  Holder’s  inequality 

{ddf  ,<t^  <c{^k—k  11^2(0)  +  11^  ~  ^  llL2(f2))  ll'/’llwo-’^(n)  ^  O’ 

as  (5  J,  0.  Thus,  is  compact  in  W~^'^{Lt),  for  each  fixed  t. 

^  J  £,(5>0 
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Multiplying  (4.2)  by  fu{k^ yields 


f{k\  +  F  {k\  +  H  {k^,  f,  =  F/  +  F/  +  F/  +  F/  +  F/ , 

where 

IV'  =  {eF/fu{k^F’^))^  +  {sufUk^F’^))^  , 

I"/  =  -e{F/f  fUk^F'^)  -e{F/f  fUk^F'^), 

=  -eF/ fuk{k^ ,F'^)ki  -  fuk{k^ ,F’^)ky, 

F/  =  Fk{k\  F'^)ki  +  Hk{k\  f,  F'^)kl  +  Hi{k\  f, 
if  =  -Uk\F’^)fk{k^F’^)ki  -  Uk^,F’^)gi{l\F’Yy 

Hence,  there  is  a  natural  decomposition  of  Cf  into  six  parts.  We  name  the  six  parts  Cf,  Ff, 

f^e,S  ^e,S  ,  ^e,5 

Regarding  Ff^ 


<C  uf 


L“(R+;L1(K^)) 


ll'^llL°°(n)  >  (/>  €(70(^2), 


which  yields  Ff  <  C,  where  A4{fl)  =  (C'c(n))*  denotes  the  space  of  bounded  measures  on  ft. 

In  view  of  Lemma  4.3  and  the  uniform  boundedness  of  the  solutions,  it  follows  that  I  F%'\  >  is 

I  ’  Je,<5>0 

compact  (and  in  fact  converges  to  zero)  in  W“^’^(H)  and  Ff  ^  ^  <  C',  for  each  fixed  t  >  0. 

Next,  for  any  (j)  G  Cc(n),  observe  that 

ff  ifcjxlxdy 
JJW 

1  _! 

l^^fl l^flj 


e  (uf)'^  dx  dy 


The  point  here  is  to  have  e  and  S  in  balance,  so  that  we  can  ensure  \sk^  \  ,  \£ky  \  <  C.  More  precisely,  we 
have  |eA:f  I  ,  \sky  \  FiC^,  and  by  choosing  e,5  according  to  (4.3)  we  achieve  this  balance.  Consequently, 


JJ <  C'll</>llLo°(n) . 


and  thus  Ff 


<  C,  for  each  fixed  t  >  0. 


Finally,  using  the  BV  regularity  of  the  coefficients  and  the  boundedness  of  the  solutions. 


J  J  ^  ‘^y  ^  C'll'/>llL“(n)  • 


and  thus  LCS’^II  ^  F,  for  each  fixed  t  >  0. 

II  ■  llA4(n) 

Similarly,  <f)dxdy  <  C  ||</>||ioom)i  thus  ^  F,  for  each  fixed  t  >  0. 

'  '  ’  A4(n) 

Summarizing,  we  have  shown  that  the  sequence  of  distributions  satisfies  the  following  two 

properties:  {i}  each  distribution  is  the  sum  of  two  terms  —  one  is  compact  in  and  the  other 
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one  is  bounded  in  In  addition,  Lemma  4.1  implies  that  {ii}  belongs  to  a  bounded 

subset  of  We  now  appeal  to  Murat  lemma  [28],  which  guarantees  that  belongs  to 

a  compact  subset  of  This  concludes  the  proof  of  the  first  part  of  the  lemma,  since  was  an 

arbitrary  bounded  open  subset  of  .  The  second  part  of  the  lemma  can  be  proved  in  a  similar  way.  □ 

Concluding  the  proof  of  Theorem  2.1.  By  combining  Lemmas  3.2  and  4.4,  we  conclude  that  •,  f)  is 

pre-compact  a.e.  for  each  t  G  [0,T].  Together  with  a  diagonal  argument,  we  can  prove  that 
converges  along  a  subsequence  a.e.  in  x  K_|_  and  in  for  each  fixed  f  >  0.  Lemma  4.2  implies 

that 

<  Ct,  Vr  G  (0,T-  r), 

and  using  this  time  continuity  estimate  it  takes  a  standard  density  argument  to  show  that  there  exists 
a  subsequence  of  that  converges  to  a  limit  function  u  a.e.  in  x  K_|_  and  in  x  K+). 

Moreover,  the  limit  u  belongs  to  L°°(R^  x  ]R_|.)  n  Lip(R+;  L^(R)). 

Equipped  with  the  strong  convergence  it  is  easy  to  prove  that  the  limit  u  is  a  weak  solution.  Moreover, 
since  each  satisfies  an  contraction  principle,  it  follows  that  the  same  is  true  for  limit.  This 
concludes  the  proof  of  Theorem  2.1.  □ 

We  close  the  paper  by  a  couple  of  remarks.  First,  we  note  that  the  apriori  bounds  in  Lemmas  4.1,  4.2 
and  4.3,  being  uniform  in  time,  enabled  us  to  deduce  pre-compactness  at  each  fix  t  >  0,  thus  circumventing 
the  temporal  argument  required  in  [36,  Appendix  A].  Second,  we  have  herein  exclusively  dealt  with 
problems  that  are  spatially  two-dimensional.  A  possible  strategy  for  going  beyond  two  dimensions  might 
be  to  adapt  the  compactness  framework  of  Panov  [30,  32,  31]. 
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